Accuracy of unperturbed motion of particles in a gyrokinetic semi-Lagrangian code

نویسندگان

  • Guillaume Latu
  • Virginie Grandgirard
  • Jérémie Abiteboul
  • Morgane Bergot
  • Nicolas Crouseilles
  • Xavier Garbet
  • Philippe Ghendrih
  • Michel Mehrenberger
  • Yanick Sarazin
  • Hocine Sellama
  • Eric Sonnendrücker
  • David Zarzoso
چکیده

Inaccurate description of the equilibrium can yield to spurious effects in gyrokinetic turbulence simulations. Also, the Vlasov solver and time integration schemes impact the conservation of physical quantities, especially in long-term simulations. Equilibrium and Vlasov solver have to be tuned in order to preserve constant states (equilibrium) and to provide good conservation property along time (mass to begin with). Several illustrative simple test cases are given to show typical spurious effects that one can observes for poor settings. We explain why Forward Semi-Lagrangian scheme bring us some benefits. Some toroidal and cylindrical GYSELA runs are shown that use FSL. Key-words: gyrokinetic, conservation laws, FSL ∗ CEA Cadarache, IRFM bat. 513, 13108 Saint-Paul-les-Durance † University of Strasbourg, 7 rue Descartes, 67000 Strasbourg & INRIA CALVI ‡ INRIA IPSO & IRMAR, 263 av Général Leclerc, 35042 Rennes ha l-0 07 27 11 8, v er si on 1 3 Se p 20 12 Précision du mouvement de particules non perturbées dans un code gyrocinétique Résumé : Une description imparfaite de l’équilibre peut conduire à des artéfacts numériques dans des simulations turbulentes gyrocinétiques. Aussi, le solveur de Vlasov and le schéma d’intégration en temps ont un impact fort sur la conservation de quantités physiques, notamment lorsque les simulations sont en temps long. L’équilibre et le solveur Vlasov doivent être finement choisis et paramétrés pour conserver les états constants (équilibre) et pour autoriser de bonnes propriétés de conservation en temps (en commençant par la conservation de la masse). Plusieurs cas tests illustratifs sont donnés pour montrer les problèmes numériques typiques que l’on peut observer si les choix pris ne sont pas adéquats. Nous expliquons pourquoi le schéma FSL (Forward Semi-Lagrangian) apporte une réponse à certains problèmes. Des simulations en configuration cylindrique et torique de GYSELA sont présentées qui utilisent FSL. Mots-clés : gyrocinétique, lois de conservation, FSL ha l-0 07 27 11 8, v er si on 1 3 Se p 20 12 Unperturbed motion of particles 3

برای دانلود متن کامل این مقاله و بیش از 32 میلیون مقاله دیگر ابتدا ثبت نام کنید

ثبت نام

اگر عضو سایت هستید لطفا وارد حساب کاربری خود شوید

منابع مشابه

A Numerical Method for Parallel Particle Motions in Gyrokinetic Vlasov Simulations∗)

A semi-Lagrangian scheme is applied for the first time to computations of charged particle motions along magnetic field lines, to numerically solve the δf gyrokinetic equations in a flux tube geometry. This new solver adopted in the gyrokinetic Vlasov simulations has an advantage over the conventional Eulerian codes in calculating the parallel dynamics, because semi-Lagrangian schemes are free ...

متن کامل

GYSELA , a full- f global gyrokinetic Semi-Lagrangian code for ITG turbulence simulations

This work addresses non-linear global gyrokinetic simulations of ion temperature gradient (ITG) driven turbulence with the GYSELA code. The particularity of GYSELA code is to use a fixed grid with a Semi-Lagrangian (SL) scheme and this for the entire distribution function. The 4D non-linear drift-kinetic version of the code already showns the interest of such a SL method which exhibits good pro...

متن کامل

LIGKA: A linear gyrokinetic code for the description of background kinetic and fast particle effects on the MHD stability in tokamaks

In a plasma with a population of super-thermal particles generated by heating or fusion processes, kinetic effects can lead to the additional destabilisation of MHD modes or even to additional energetic particle modes. In order to describe these modes, a new linear gyrokinetic MHD code has been developed and tested, LIGKA (Linear Gyrokinetic Shear Alfvén Physics) [1; 2], based on a gyrokinetic ...

متن کامل

A 5D gyrokinetic full-f global semi-Lagrangian code for flux-driven ion turbulence simulations

This paper addresses non-linear gyrokinetic simulations of ion temperature gradient (ITG) turbulence in tokamak plasmas. The electrostatic Gysela code is one of the few international 5D gyrokinetic codes able to perform global, full-f and flux-driven simulations. Its has also the numerical originality of being based on a semi-Lagrangian (SL) method. This reference paper for the Gysela code pres...

متن کامل

Improving conservation properties of a 5D gyrokinetic semi-Lagrangian code

In gyrokinetic turbulent simulations, the knowledge of some stationary states can help reducing numerical artifacts. Considering long-term simulations, the qualities of the Vlasov solver and of the radial boundary conditions have an impact on the conservation properties. In order to improve mass and energy conservation mainly, the following methods are investigated: fix the radial boundary cond...

متن کامل

ذخیره در منابع من


  با ذخیره ی این منبع در منابع من، دسترسی به آن را برای استفاده های بعدی آسان تر کنید

برای دانلود متن کامل این مقاله و بیش از 32 میلیون مقاله دیگر ابتدا ثبت نام کنید

ثبت نام

اگر عضو سایت هستید لطفا وارد حساب کاربری خود شوید

عنوان ژورنال:
  • CoRR

دوره abs/1209.0317  شماره 

صفحات  -

تاریخ انتشار 2012